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Abstract 



We describe the quantitative refinement of nanopar- 
ticle structures from gold nanoparticles probed by ul- 
trafast electron crystallography (UEC). We estabhsh 
the equivalence between the modified radial distribu- 
tion function employed in UEC and the atomic pair 
distribution function (PDF) used in x-ray and neu- 
tron powder diffraction analysis. By leveraging PDF 
refinement techniques, we demonstrate that UEC 
data are of sufficient quality to differentiate between 
cuboctahedral, decahedral and icosahedral nanopar- 
ticle models. Furthermore, we identify the signatures 
of systematic errors that may occur during data re- 
duction and show that atomic positions refined from 
UEC are robust to these errors. This work serves as a 
foundation for reliable quantitative structural analy- 
sis of timc-resolvcd laser-excited nanoparticle states. 



1 Introduction 

Ultrafast electron crystallography (UEC) is a promis- 
ing new approach for studying the structure of 
photoexcited materials on pico and femtosecond 
timescales.[l] [2] It extends existing ultrafast pump- 
probe spectroscopy techniques by providing vital 
structural information on the photoexcited states, 
which is required to obtain a more complete under- 
standing of these transient states. It has long been a 
dream to make "molecular movies" [3l |4j [5] of atoms 
as a material transforms from one state to another. 
The use of x-rays for this purpose, obtained from 
powerful synchrotron sources and in the future from 
hard x-ray free electron lasers, has the advantage that 
the data can be modeled quantitatively using existing 
x-ray crystallographicfB] and atomic pair distribution 
function]?] (PDF) approaches. However, on the ex- 
perimental side, the use of electrons as a probe has 
a number of advantages over synchrotron x-ray tech- 
niques. First, the scattering power of the electrons is 
higher, allowing smaller and more dilute samples to 
be studied. Second, it is easier to control the pulse 
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structure of the probe-pulse on femtosecond time- 
scales without the need for fast choppers. [8] 

A challenge when using electrons as a probe is that 
the scattering is strong, resulting in increased mul- 
tiple scattering, and the single-scattering kinemati- 
cal approximation used in most x-ray analyses may 
not be valid. Recent developments in quantitative 
electron crystallography, [S] especially with the devel- 
opment of precession methods of data collection. [TU] 
show that these problems can be overcome under cer- 
tain circumstances and quantitative structure solu- 
tions are possible using electrons. Another challenge 
is structure solution when the underlying structural 
correlations only extend over nanometer dimensions, 
the so-called "nanostructure problem" . [TT] Here we 
explore whether methods developed for quantitative 
structure solutionjT^] and refinement [T^ UM ITS] of x- 
ray or neutron derived PDFs can be extended to UEC 
data from nanoparticles. The full potential of UEC 
can be realized only when the quantitative reliability 
of structures determined from the data is established. 
We investigate this issue here for the case of data from 
metallic nanoparticles tethered to a surface. 

There is a large literature on the structure solu- 
tion of small molecules in the gas phase from electron 
diffraction data.[16j However, the UEC experiments 
of interest are on samples in the solid or liquid phase 
where the validity of the kincmatical approximation 
is less clear. Quantitative structure solutions from 
electron crystallography based on kinematic scat- 
tering are possible using very thin crystals, [T71 fT5] 
weakly scattering specimens |19l 120] and nanoparti- 
cles containing just 10-20 unit cells. [H] In these cases, 
the probability of a multiple-scattering event is suffi- 
ciently small that it can be treated as a minor pertur- 
bation on the single-scattering signal, and the data 
can be properly corrected for a quantitative struc- 
ture solution. Here we investigate whether UEC 
data from nanoparticles tethered to a surface can be 
treated in this way. There are powerful tools emerg- 
ing for the study of nano-scale structures from ensem- 
bles of nanoparticles using x-ray and neutron diffrac- 
tion data, [ni [H m m [55] These tools can also 
be applied and extended for studying UEC data if the 
quantitative reliability of the data in the kinematical 
limit can be established. 



In this paper we show that it is possible to ob- 
tain UEC data of sufficient quality for quantita- 
tive structure modeling. We investigate the robust- 
ness of structural solutions with respect to uncon- 
trolled aspects of the data reduction such as the use 
of empirical methods for removing incoherent back- 
grounds and determining an effective electron-form- 
factor that considers the absorption by the substrate. 
We elucidate potential sources of error in the data re- 
duction procedure, and indicate the signature those 
errors. Additionally, we show the relationship be- 
tween the correlation functions traditionally used in 
gas-phase electron diffraction, and those used in the 
study of nanostructure with x-rays and neutrons, 
bringing together these two fields of study. These re- 
sults place the UEC method on a more sound footing 
with respect to the study of ensembles of nanoparti- 
cles, opening the door to reliable quantitative struc- 
ture solution of transient states, and therefore physi- 
cally meaningful molecular movies of these materials. 

2 UEC with nanoparticles 

The photochemical processes of bond breaking and 
forming can be resolved using femtosecond laser tech- 
niques in a pump-probe configuration, in which the 
first laser pulse (pump) is used as a trigger to ini- 
tiate the chemical reactions and the ensuing probe 
laser pulse serves to monitor the changes via spectro- 
scopic responses. [3] By replacing the probe laser pulse 
with a laser-triggered ultrashort electron pulse, the 
transient structural dynamics can be directly probed. 
This ultrafast electron diffraction (UED) [27] method 
was shown to be especially advantageous in the stud- 
ies of optically dark processes, such as nonequilib- 
rium structural transformation. [28] and radiationless 
relaxation channels for the photo-excited state. [29] 

The UEC development extends UED to the con- 
densed phase, and has been applied to study the 
collected excitations involving photons, electrons 
and phonons from the interfacial water bilayers 
to superconductors . [501 [5T1 [32] Whereas the peri- 
odic orders from the crystalline samples offer higher 
signal-to- noise ratio (SNR), and potentially, higher 
spatiotemporal resolutions, the complexity of the 
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anisotropic scattering in the dynamical scattering 
regime poses significant challenges in quantitatively 
assessing all the possible details of the relevant atomic 
dynamics that the higher SNR can now provide. 

Recently, UEC has been applied to investigate 
nanoparticles that are highly dispersed on a sur- 
face to produce powder-like patterns for size-scale 
down to 20 A. [2] To suppress the background signals 
from the supporting substrate, soft-anchoring of the 
nanoparticles using self-assembled molecular (SAM) 
layer on silicon surface is employed. The implemen- 
tation of this specific UEC for studying nanoparticles 
(as shown in Fig. [T|) affords a low coverage to ensure 
the scattering from isolated particles can be recorded 
unobstructively, and provides the needed chemical 
and physical isolation of nanoparticles under repeti- 
tive time- resolved studies. [33] 

Even though a relatively clean coherent diffrac- 
tion pattern for nanoparticles can be obtained this 
way by rendering the scattering from the substrate 
largely diffusive, it is still unclear how the residual ef- 
fects caused by the uncontrolled diffusive background 
and the surface absorption affect the quantitative re- 
trieval of structural functions. The empirically de- 
rived structural functions from UEC|33^ will be re- 
fined using standard x-ray PDF analysis to examine 
the self-consistency and robustness of UEC data for 
quantitative structural analysis. 



3 Correlation functions 

In this section we discuss the relationship between 
the modified radial distribution function (mRDF) 
obtained from UEC experiments with the PDF ob- 
tained from x-ray and neutron measurements. His- 
torically, the mRDF has been applied in gas-phase 
electron diffraction to study atom-pair correlations 
in small molecules containing a few to tens of atoms. 
In UEC, the mRDF is used to study the change in 
atomic arrangements in excited-state bulk- and nano- 
structures.pilM] 

The mRDF, R{r), is defined as the trun- 
cated Fourier transform of the molecular scattering 




Figure 1: UEC for nanoparticles experimental setup, 
implementing the laser-pump-electron-probe config- 
uration. The electron pulse is delayed relative to the 
'pumping' laser pulse in a time-sequence (At) to ex- 
amine the structural dynamics. The incident angle 
of the electron beam is kept low, typically 1-2°, to 
achieve high sensitivity to the nanoparticles, which 
are dispersed on a self-assembled molecular mono- 
layer atop a supporting silicon substrate, to produce 
powder-like diffraction patterns. 
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intensity. [55] 

R{r) ^ sM{s)sm{sr) ds, (1) 

Jq 

where s is the magnitude of the scattering momentum 
transfer and Smax is the maximum determined value 
of s. The molecular scattering intensity, M{s), con- 
tains the structural information extracted from the 
scattered electron intensity. The mRDF and sM{s) 
contain the same structural information, and it is 
the sM(s) curve that is typically used in structural 
modeling. [TB] 

The PDF, G{r), is used in x-ray and neutron inves- 
tigations of liquids, and amorphous and nanostruc- 
tured materials. [7] The PDF is the truncated Fourier 
transform of the reduced total scattering structure 
function, F(Q).[37] 

G{r) ^ - r^^'^ F{Q) sin(Qr) dQ, (2) 

where Q is the magnitude of the scattering mo- 
mentum transfer [38] and Qmin and Qmax are the 
measured Q-extrema. For elastic scattering, Q = 
(47r/A) sin(6'), where A is the wavelength of the probe, 
and 26 is the scattering angle. Note, the scattering 
angle in the gas phase electron diffraction literature is 
generally denoted 0, rather than 26. Notation aside, 
Q and s are equivalent: 

Q = s. (3) 

In deference to the respective literatures, and to help 
differentiate between the two sets of functions, we 
will use s when referring to UEC and Q in reference 
to x-ray and neutron diffraction. 

As we will show in the next section, F{Q) and 
sM{s) are effectively the same, where the only dif- 
ferences come about due to experimental effects and 
uncertainties in the data reduction process. These 
being equal, the important difference between the 
mRDF and PDF is the lower integration limit in their 
respective definitions. |37| 

There is always a finite lower bound on the mea- 
sured scattering momentum (s or Q) due to experi- 
ment geometry and physical limitations of measure- 



ment equipment. It is the practice in gas-phase elec- 
tron diffraction and UEC to compensate for this miss- 
ing scattering intensity by using a presumed struc- 
tural model.[35l [2] In gas-phase electron diffraction 
the missing intensity below Smin is extrapolated to 
s = using the model, and in UEC the same effect 
is achieved by adding a baseline to the transformed 
real-space signal. This works well when the system 
under study is made up of individual small molecules 
in a dilute gas phase. However, it becomes highly 
problematic in condensed systems. In this case, the 
convention in x-ray and neutron diffraction is to ig- 
nore any scattering in the low-Q limit. 

When the Fourier transform includes the intensity 
extending to s = Q = 0, one obtains the mRDF, 
which is always positive (apart from termination ef- 
fects due to a finite Smax)- However, if the small 
angle scattering is excluded, one obtains the PDF, 
which is a function that oscillates around zero, sit- 
ting on top of a negative baseline. [37] In either case, 
with the same minimum and maximum momentum 
transfer, the same structural information can be ob- 
tained from the mRDF or the PDF. We use the PDF 
in this study in order to avoid the extra step of com- 
pensating for the missing scattering information. 

4 Data reduction in detail 

In this section we will describe how sM[s) is obtained 
in practice in order to examine the systematic errors 
that may affect the extraction of structural informa- 
tion from the data. We will first discuss how the total 
scattering structure function is typically obtained in 
x-ray and neutron diffraction. This process has an es- 
tablished theoretical foundation, [7] and will serve as 
a basis of comparison with the UEC data reduction 
procedure. 

4.1 X-ray and neutron diffraction 
data reduction 

The reduced total scattering structure function, 
F{Q), is defined in terms of the total scattering 
structure function, S{Q), as F{Q) = Q{S{Q) - 1). 
The structure function contains the discrete coher- 
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ent singly scattered information available in the raw 
diffraction intensity data. It is defined according 
to[37j 

iciQ) {{.f-im 



SiQ) 



(4) t 



which gives 



5(g) - 1 



IciQ) - N{p 
idiQ) 



(5) 



where / is the Q-dependent x-ray scattering factor 
or Q-independent neutron scattering length as appro- 
priate and (. . .) represents an average over all atoms 
in the sample. In this equation, IdQ) is the coher- 
ent single-scattered intensity per atom and Id{Q) is 
the discrete coherent scattering intensity, which ex- 
cludes the self-scattering, iV(/^).[37j The coherent 
scattering intensity is obtained from the measured 
intensity by removing parasitic scattering (e.g., from 
sample environments), incoherent and multiple scat- 
tering contributions, and correcting for experimen- 
tal effects such as absorption, detector efficiencies, 
detector dead-time and so on. [7] The resulting cor- 
rected measured intensity is normalized by the inci- 
dent flux to obtain IdQ)- The self-scattering, N{p), 
and normalization, N{f)^, terms are calculated from 
the known composition of the sample using tabulated 
values of /. 

As evident in Eq. [H to obtain S{Q) — 1 from 
IciQ) "we subtract the self-scattering, N{f^), which 
has no atom-pair correlation information, and divide 
by iV(/)^. As a result, S{Q) — 1 oscillates around 
zero, and asymptotically approaches it at high Q as 
the coherence of the scattering is lost. If the exper- 
imental effects arc removed correctly, the resulting 
F{Q) and G{r) are directly related to, and can be cal- 
culated from, structural models. |2Z] The corrections 
are well controlled in most cases and refinements of 
structural models result in reduced values that 
approach unity in the best cases. Some uncertainty 
in the corrections can be tolerated because they are 
mostly long-wavelength in nature (for example, the 
absorption correction), which results in aberrations 
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Figure 2: The UEC data reduction procedure show- 
ing the reduction of I{s) to M{s). 

in G(r) at very low-r below any physically meaning- 
ful region of the PDF. [39] 

4.2 UEC data reduction 

In practice, M{s) is obtained from the raw experi- 
mental intensity, /(s), by an empirical approach. |33| 
This approach is necessary because the reflection ge- 
ometry employed to effectively sample the surface- 
dispersed nanoparticles results in a shadowing effect, 
i.e. the diffraction intensities collected at small exit 
angles are strongly attenuated due to surface absorp- 
tion. In practice, due to the random nature of the sur- 
face roughness, the attenuation factor is not known 
precisely. It is not immediately clear that the empiri- 
cal method results in quantitatively reliable mRDFs, 
as demonstrated in a more simple transmission ge- 
ometry in which the tabulated atomic scattering form 
factor can be used directly to yield quantitative struc- 
ture solutions of molecular systems. [23 

The UEC data reduction procedure is demon- 
strated in Fig. [21 First, an empirical background is 
estimated by fitting a slowly varying function, such 
as a low order polynomial, through the center or base 
of the raw intensity. This background is subtracted 
from the raw intensity. The background-subtracted 
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data are then divided by the same background func- 
tion. We denote this estimated background as Ib{s) 
and can write the resulting intermediate intensity, 
M'{s). as 



M'{s) 



his) 



(6) 



The next step in processing the data is to adjust for 
any offset in the initial background. This secondary 
background, M^(s), captures the slowly- varying com- 
ponents of M'(s), and as a result the remaining sig- 
nal oscillates around zero. The result is the molecular 
scattering intensity, 



M{s) = M'{s) - Ml{s). 



(7) 



4.3 Practical relationship between the 
molecular scattering and the total 
scattering structure function 

In order to understand how M{s) relates to S{Q), we 
consider an ideal case; scattering from a monatomic 
system with no undesirable experimental effects in- 
fluencing the data. If measured without an offset. If, 
captures only the slowly varying self-scattering from 
the (normalized coherent) scattering intensity. This 
self scattering is equal to N{f)'^ for a monatomic 
system. Thus, M'{s) = - N{p))/N{ff = 

Id{s)/N{f)^ ~ S{s) — 1. Since this background is 
measured perfectly, A/^ = 0, and so M{s) = M'{s) = 
Sis) - 1. 

In reality, the data are not ideal; must correct for 
experimental effects, and Af^ must somehow compen- 
sate for the accuracy of More formally, we write 
the raw intensity as the coherent intensity modified 
by multiplicative term that captures multiplicative 
experimental effects (e.g., absorption correction) and 
an additive term representing the incoherent back- 
ground and additive experimental effects (e.g., sub- 
strate scattering, inelastic scattering) 



lis) =a(s)/c(s)+i^(s). 



(8) 



The background intensity, his), estimates the 
modified self-scattering (t(s)7V(/^), the incoherent 
scattering background, and the undesirable additive 
effects, i^is). The goal in subtracting his) is to 



remove the self-scattering and additive effects from 
the intensity. This background is typically estimated 
using a low-order polynomial, which can approxi- 
mate the experimentally modified iV(/^) with rela- 
tive accuracy. [30] The estimate is subject to error, 
however, which may be due to components of i'(s) 
that are poorly estimated or due to an intentional 
offset of the estimated background curve. We denote 
this error by — e(s). Thus, the first stage of the UEC 
data reduction results in 



M'is) ^ 



ais)his) + ^(.s) - iais)N{p) + i^js) ^ ejs))) 
(7(s)7V(/2) + i.(s) - e(s) 
<jis)ihis)-N{p)) 



+ 



a(s)iV(/2) + i.(s) - e(s) 
eis) 

ais)N(p)+:,{s)~eis) 

<jis)ihis)-N{p)) 
a(s)iV(/2) + i.(s) - e(s) 



lis)- 



(9) 



Note that 7(5) and the denominator in Eq. [2] are 
slowly varying functions of s. 

The next step is to eliminate the effects of 7(5), 
which result from the offset in the estimate of his)- 
Since 7(3) is slowly varying it can be fit with a low- 
order polynomial. With 7(5) ehminated, the result 
will oscillate around zero. There will be some error in 
this estimate as well, which we denote ~(3{s). Thus, 



Miis) = 7(s) - Pis). 
Using this expression for Ajf^(s) gives 



(10) 



ais)ihis)-N{p)) , ^ 
a[s)N{f^} + v{s) - e(s) 

This is very near to the form of Eq. [5l To get this in 
the desired form, wc write Ids) — N{p) ~ (5'(s) — 
l)iV(/)2 (sec Eq.© so that 



Mis)=ais)iSis)-l)+Pis), 



where 



a(s) 



ais)N{fy^ 



ais)N{p)+uis)-eis)' 



(12) 



(13) 
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and we note that a{s) is also slowly varying. Now, 
understanding the effective difference between M (s) 
and S{s) — 1 reduces to understanding the effects a{s) 
and (3{s) have on the measured correlation function. 

First we consider (3{s). By construction, (3{s) is 
small, and contains only slowly varying components. 
As mentioned above, this means that (3{s) contributes 
only to the low-r region of the correlation function, 
below the physically interesting region. [39] If the ex- 
perimental effects are not slowly varying, then the 
Fourier transform of /3(s) will leak into the struc- 
turally relevant portion of the correlation function. 
The effect is unpredictable, but should be small. 

The multiplicative term, a(s), is likely to have a 
more significant effect on the correlation function. In 
the best case, a{s) is constant and it scales the peaks 
of the correlation function uniformly, which does not 
obscure the structural information. The only aspect 
of a{s) that is controlled by the data reduction is 
how quickly it oscillates. This is controled by how 
well the modified self-scattering is fit by the inten- 
sity background (manifest in i/{s) and e(s)). De- 
pending on how well these factors are estimated, a{s) 
may destort the peak profile of the correlation func- 
tion. At worst, this could result in a measurable peak 
shift, which would complicate the extraction of reli- 
able structure parameters. 

5 Methods 

Data have been collected on size-selected gold 
nanoparticles with diameter of approximately 
20 A using the UEC experimental setup described 
in ^ Further details on the sample preparation 
and data collection can be found elsewhere. To 
gain quantitative information about the morphology 
of gold nanoparticles we modeled PDFs from ex- 
perimental UEC data using three structure models: 
a 309-atom cuboctahedron (diameter 22.3 A), a 
309-atom Mackay icosahedron[^ (diameter 23.1 A) 
and a 181-atom decahedron [32] (diameter 23.8 A). 
These gold cluster types have been shown to exist at 
a similar size|43| and have been the subject of a pre- 
vious UEC study. [33] We also fit these nanoparticle 
models with one fewer and one additional shell. The 
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Figure 3: PDFs from likely models for ~ 20 A gold 
nanoparticles. Top: Cuboctahedral, decahedral and 
icosahedral models. Bottom: PDFs calculated from 
the cuboctahedral (top), decahedral (middle) and 
icosahedral (bottom) models. 

results in those fits were considerably worse than for 
these nanoparticle sizes, and will not be reported. 

The structure models are shown in FigJ3] along 
with the theoretical PDF from each modellll As can 
be seen in that figure, the models give PDFs that 
are very different past 5 A, which reflects their topo- 
graphical differences. All nanoparticles were modeled 
with rigid atomic arrangements and uniform isotropic 
atomic displacement parameters (ADPs). 

In all reflnements the nanoparticle models were al- 
lowed to expand or contract isotropically keeping the 
number and relative positions of atoms in each model 
fixed. The samples are nominally mono-disperse, so 
to avoid complicating the analysis we have not at- 
tempted to model size distributions. Furthermore, 



^Figures created with the AtomEye structure viewer. [44] 
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we have not considered strain effects besides a simple 
isotropic expansion of the models. 

We calculate F{Q) from a structure model us- 
ing the Debye equation, [ISl [13] modified such that 
the pair-contributions are attenuated by a Dcbye- 
Waller factor, and this is transformed to G(r) ac- 
cording to Eq. [H The Debye- Waller factor is derived 
from the ADPs and a refinable vibrational correla- 
tion term. [JSl |37] A scale factor was refined to ac- 
count for the undetermined scale of the raw intensity. 
Resolution factors |3B] that broaden and dampen the 
PDF peaks were also refined to simulate the effects 
of the finite s-resolution of the measurement. The re- 
finements were performed using a home-written least 
squares regression algorithm based on the PDFfit2 
program. [14j 

We distinguish the fits by their agreement with the 
data. This is quantified with the weighted residual of 
the PDF, defined as 



/ E.^(r-.)(G(r.)-G^(r,))^ 
E,c.(rOG(r,)2 



(14) 



where G denotes the experimental PDF and Gm is 
the PDF calculated from a structure model. For this 
study the weights tj(ri) are equal for every point and 
therefore cancel in R-u, . The factor was calculated 
only over the physical portion of the fit range, which 
we estimate to be between 2.2 and 20 A. 

To test the robustness of the data reduction and 
refinement procedure, we have also performed re- 
finements on the PDF produced from four different 
data reduction protocols and compared the results. 
These protocols differ in how the intensity back- 
ground, /f,(s), is determined. The different protocols 
produce different s-dependent scaling factors, a{s), 
and ultimately different PDFs. For two reduction 
procedures. Ib{s) is determined by fitting along the 
bottom of I{s) using 4th and 5th degree polynomials 
(denoted poly4b and polyBb , respectively), the lat- 
ter of which is the standard procedure. [T] The other 
two procedures fit Ib{s) through the center of I{s) us- 
ing a 7th degree polynomial (denoted poly/c ) and an 
exponential function of the form y^s — sq exp(— s'^'''^) 
(denoted expn) that is empirically found to well de- 
scribe the shadowing effect and decay of the scatter- 
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Figure 4: Fits to the standard data. polySb data 
(blue dots), fit (red line) and difference (green line, 
offset). (Top) Cuboctahedral model. (Middle) Deca- 
hedral model. (Bottom) Icosahedral model. Note the 
comparative quality of the cuboctahedral fit above 
5 A and generally poor fit to the third and fourth 
peaks. 



ing signal. 



6 Results 

6.1 Quantitative Structure Modeling 

To determine whether the UEC data are of sufficient 
quality to differentiate, quantitatively, between struc- 
tural models, we fit cuboctahedral, decahedral and 
icosahedral nanoparticle models to UEC PDF data 
produced using the standard procedure {poly5b ) as 
described in ^ 

Among the models used in these refinements, the 
cuboctahedral model stands out as the best, result- 
ing in a Rw value of 0.23, compared to 0.30 for the 
decahedral model and 0.41 for the icosahedral model. 
These fits fits are shown in Fig. ID None of the struc- 
ture models reproduce the third and fourth peaks of 
the PDF well. It is possible that this misfit is due 
to errors introduced during the data reduction. We 
explore this possibility below. The structural results 
from the cuboctahedral model are shown in Table [T] 
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6.2 Robustness 



Given the empirical nature of the methods for pro- 
ducing the PDF from the diffraction intensity, it is 
worthwhile to investigate the robustness of refine- 
ment results based on the different procedures. This 
will give insight into the influence of the data re- 
duction procedure on the PDF, and will give us a 
quantitative measure of the systematic bias that may 
be introduced during data reduction. We have per- 
formed refinements of the cuboctahedral model de- 
scribed above to the PDF produced from the four 
different data reduction procedures described in fJS] 
and demonstrated in Fig. [S] 

As can be seen in Fig. O the data reduction proce- 
dures all produce sM{s) and PDF curves similar to 
each other. The major features of these curves (indi- 
cated by the dashed lines in the figure) are preserved 
using each procedure. The largest differences among 
the PDFs appear in the peak widths and amplitude 
ratios. Peak asymmetry is hard to detect given that 
the PDF sits on a negative baseline, and the peaks 
are therefore not Gaussian. The bottom panel of 
Fig. [S] shows the difference between the PDF from 
the polySb background (the standard procedure) and 
the other PDFs, scaled so that the heights of the first 
peaks are identical. The pre-peak signal is rather dif- 
ferent in each signal, indicative of the additive errors, 
/3(s), in sM{s) resulting from the data reduction pro- 
cedures. 

Fits to the data were performed using the cubocta- 
hedral model as described in SjSl The cuboctahedral 
model fits can be seen in Fig. [5] and are summarized 
in Table[TJ In the table, a is the effective FCC lattice 
constant derived from the interior nearest neighbor 
distance. As seen from the figure and table, all fits 
are of variable quality, with a markedly improved Ru, 
value for the exponential background. It can be seen 
that this fit reproduces the fourth peak significantly 
better than the other fits. 

From Table [1] we see that the measured peak width 
parameter (Uiso) is not robust with respect to the 
data reduction procedure. This indicates that the at 
least some of the data suffer from peak profile dis- 
tortions due to the data reduction protocol. Con- 
sidering the ideal data reduction scenerio discussed 




r(A) 

Figure 5: Similarity between structure functions 
from various data reduction procedures. Blue line: 
poly4b . Red line: poly5b . Green line: polyYc . 
Brown line: expn. See text. Note the similarity of 
the major profile features, as indicated by the dashed 
lines, (a) Ib{s) for each protocol through the raw in- 
tensity (black circles), (b) sM{s) for different re- 
ductions, (c) PDFs from the different reductions, 
(d) PDF differences taken from the polySb PDF. AU 
PDFs are scaled to have the same first-peak ampli- 
tude. The structural region is r > 2.2 A. 
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Table 1: Refinement parameters from cuboctahe- 
dral fits. Here, a is the effective lattice parameter 
derived from nearest neighbor distances and Uiso is 
the isotropic ADP (see text) . No standard deviations 
are shown on the refined parameters as we don't know 
uncertainties reliably on the raw intensities. 



in ^4.31 where the ideal intensity background per- 
fectly matches the modified self-scattering through 
the center of the raw intensity, we believe that proto- 
cols poly/c and expn, which estimate the background 
through the center of the raw intensity, introduce less 
error into the correlation function. This conclusion 
is supported by the superior low-r fit agreement to 
these data and the smaller fit residuals {Rw)- 

The effective lattice constant (a) from Table [T] is 
robust with respect to the data reduction procedure. 
This shows that the errors in the correlation function 
introduced by the data reduction procedures are not 
severe enough to result in a measurable peak shift. 
Thus, the geometric structure information is read- 
ily accessible from the UEC data despite the uncon- 
trolled aspects of the data reduction protocol. 

We note that the estimated lattice parameter is 
smaller than that for bulk gold. This is not a physi- 
cal effect, or a result of the data reduction protocol, 
but rather due to a lack of instrument calibration. 
Driven primarily by the interests of differential struc- 
tural changes, UEC frequently treats the data using 
the known ground state structure as the reference. 
In the case where the ground state structure is not 
known, only relative changes can be obtained. Given 
the very small uncertainty of the refined lattice pa- 
rameters shown here, these ground state data could 
be used to calibrate the UEC setup by comparing to 
data from the same particles measured with x-rays in 
a conventional PDF measurement, thus making the 
absolute structural solution from UEC possible. 



I ' I 
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0.0250 
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0.20 
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0.0167 


expn 


0.20 


4.116 


0.0185 
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15 



20 



Figure 6: Fits using the cuboctahedral model to 
PDFs produced with different data reduction proto- 
cols. (See Fig.[5]for color key.) PDF fits (solid lines), 
data (gray circles) and difference between data and 
fit (solid lines, offset below). 

7 Conclusions 

We have shown the equivalence of the mRDF com- 
monly used in UEC and the PDF typically employed 
in x-ray and neutron powder diffraction. This has 
allowed us to use powder diffraction modeling tech- 
niques to refine cuboctahedral, icosahedral and dec- 
ahedral models to the PDF derived from UEC data. 
We have found that the cuboctahedral model agrees 
best with the data. This demonstrates that UEC 
data can be of sufficient quality to differentiate be- 
tween different nanoparticle structure models. 

We have tested the cuboctahedral model with mul- 
tiple PDFs produced using different data reduction 
protocols that differ in how the intensity background 
is removed and the data normalized. Given the large 
differences in the background curves fit through the 
raw intensity data, we are able to determine the ef- 
fective lattice parameter of the cuboctahedral model 
with remarkable certainty. We cannot determine with 
certainty the thermal displacement parameter within 
the model because this parameter is strongly infiu- 
enced by the form of the intensity background. 

The figure of merit of our best fits (cuboctahe- 
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dral polyVc and expn, = 0.20) is comparable 
to results from high-quality nanoparticle x-ray PDF 
studies. [551 ES] The cuboctahcdral model affords lit- 
tle flexibility to compensate for systematic errors, so 
the results indicate that these errors are not signifi- 
cantly worse than what one would expect from syn- 
chrotron x-ray data. We conclude that, with careful 
data reduction and a standard system for calibration, 
UEC is applicable as a quantitative nanostructure 
probe. 
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